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ABSTRACT 

It is shown that quantum fluctuations due to a nontrivial gravitational background 
in the flat radiation dominated universe can play an important cosmological role gen- 
erating nonvanishing cosmological global charge, e.g. baryon number, asymmetry. The 
explicit form of the fluctuations at vacuum and at finite temperature is given. Impli- 
cations for particle physics are discussed. 
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1. Introduction 



Since the early eighties it has been widely recognized that quantum fluctuations of 
scalar matter fields may play an important role in cosmology, especially in the context 
of the inflationary de Sitter epoch [1]. Actually, the case of the scalar field in the de 
Sitter space has been the most extensively and carefully studied one. The reason is 
that in the de Sitter space fluctuations of the light fields (m 2 /H] « 1, where if/ is the 
de Sitter Hubble parameter) grow linearly with time assuming finally a significantly 

H 4 

large value of the order of which singles out the de Sitter universe. It is believed 
that fluctuations produced at the time of inflation are seen during subsequent stages 
of the evolution of the universe as energy density inhomogeneities responsible for the 
formation of the large scale structure. It is also advocated that those fluctuations set 
initial conditions for the classical evolution of fields in subsequent epochs. 

In contrast to the above, it is usually assumed that gravitationally induced scalar 
field fluctuations in spatially flat radiation dominated (RD) and matter dominated 
(MD) epochs are irrelevant for particle cosmology. We would like to point out that 
this assumption is not properly discussed in the literature. On one hand, one observes 
that in the RD universe the fluctuations (as explained in this letter) decrease in time. 
On the other hand, they may in principle be large enough to control violation of some 
symmetries or to alter the evolution of some fields present in field-theoretical models. 
We want to stress that this problem becomes particularly important in view of the 
ongoing search for a reliable mechanism for production of the baryon asymmetry in 
the Universe, the need of the better understanding of the scenarios for late-time phase 
transitions and discussions of the possible lepton number nonconservation. 

In this letter we would like to address explicitly the problem of quantum fluctuations 
of the massive scalar field during the RD epoch. This epoch covers most of the history 
of the universe, and the temperature range from, say, 10 14 GeV down to 10 eV. On 
that energy scale one can find a lot of interesting phenomena in popular extensions 
of the standard model such as its supersymmetric version or string inspired models, 
what, in our opinion, justifies the research reported in this work. 

The paper is organized as follows. In Section 1 we set our notation and subsequently 
evaluate fluctuations of a massive scalar field in the RD flat Robertson- Walker space 
at vacuum and at finite temperature. In Section 2 we apply our formulae to a general 
field theoretical model with particular attention paid to two specific examples resem- 
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bling the Affleck-Dine model [2] and the so called spontaneous baryogenesis scenario 
[3]. Finally, in the last Section we review our results and present conclusions. 



1. Scalar field fluctuations in the RD universe 



The RD Universe is the solution to the Einstein's equations with the energy- 
momentum tensor in the form = diag(—p,p,p,p). Tracelessness of the implies 
equation of state for the content of the RD Universe: p = 3p. In this letter we assume 
a flat RD space endowed with the metric g^ u = diag(—l,a 2 (t),a 2 (t),a 2 (t)) where a(t) 
is the RW scale factor given by a(t) = (w) 1 / 2 with u defined as u = t/t , t being the 
beginning of the RD epoch. 

We couple a massive scalar field to gravity in the minimal way ( note that here the 
curvature scalar R vanishes identically ) 

S[<j>] = J d A x^g (-g^d^d^ - m 2 2 ) . (l) 

As usual in this type of analysis we assume that there is no "back-reaction" of the 
scalar field on the metric, cf. [4]. The equation resulting from eq.(|l]) is 

where <fik( u ) is the spatial Fourier transform of the field (f>(x,t), 

m t) = J {MuyM + h.c.) (3) 

with aJ and its hermitian conjugate denoting standard creation and anihilation opera- 
tors respectively. 

A general solution of the equation (|2]) is given by confluent hypergeometric functions 

(j, k {u) = AAk, m)2ikt e- imuto 1 F 1 3/4 + °-, 3/2, 2imt u 

\ 2m J 

+ A 2 (k,m)^e- imut \F 1 ^l/A + ^,l/2,2imt uj. (4) 

The two undetermined coefficients A% and A 2 , which may in principle depend on both 
m and k, are not independent if one takes into account quantization condition imposed 
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on a field 



{x,t),d t cf>(y,t)} = -j=8W{x-y) (5) 
= (2ic) 3 2k5 ( ' 3 \k - k') (6) 



From above equations and decomposition (|3p we get a normalization condition 

lm(Mu)d t Mu)) = -j= (7) 
which translates into the constraint on A\ and A 2 ( we confine them to be real ) 

A 1 {k,m)A 2 {k,m) = 1 (8) 

In most general case there are several ways of fixing both coefficients. One possibility 
is to use initial conditions set at the timelike surface t = t for <fi and d t <p. This is 
the proper procedure if one knows for example the explicit solution for <p in the epoch 
preceding the RD one. We do not assume such a detailed knowledge, hence we use 
an alternative approach instead. We demand that the "correct" mode functions we 
choose, which will define our Fock space, should approach at short distances (k — > oo) 
the massless positive frequency solution, 

<f> k (u) -> -Le 2lkt °^ (9) 

u I 

In this way we obtain the asymptotic behaviour of both coefficients 

A 1<2 (k,m)^l (10) 

Here we assume that A\ — A% — 1, what completes the definition of our Fock space. 

We set out to calculate the fluctuations of the field <fi i.e. < O|0 2 |O >. This quantity 
is badly divergent and needs renormalization. We define the renormalized fluctuations 
as the difference between RD and the Minkowski space fluctuations, hence the relevant 
object to look at is the difference 

< 4> 2 > R =< 4> 2 > - < 4> 2 > M (11) 

As we shall see, this definition gives the finite result. It may easily be checked that 
in terms of Fourier modes 0fc(w) the renormalized fluctuations (|ll]) are given by the 
formula 
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Using (||) we can write down the explicit formula 



< C >R 



Ant* 



POO 




/ dy 


y 


lo 





2mt 2' 



\ ^ /3 n/ 2 3 

2,mfl+2.y 1 f\» i + — .-.fcmf 



y/y' 2 + m 2 t' 2 



(13) 



where y = kt/ y/u. Unfortunately, the above expression cannot be evaluated in its most 
general form. However, it is possible to write down the systematic expansion of the 
mode functions and the integral (|T2"D in terms of mt. Using such an expansion we 
will be able to discuss reliably fluctuations in the regime of small mass and to control 
the passage to the massless limit. For the region mt > 1 we will have to rely on the 
numerical calculations. 

In the case mt < 1, the relevant expansion of modes is given by (cf. [5]) 

1 00 1 
-Y. \pt 1/2) (2tmt) Jn ^(2y)+tp^ 2 H2tmt) Jn (2y)} — 
u 



n=Q 



(Zy) 71 - 1 



(14) 



M( z ) w n = e «/2(coth(2w)-l/2to) 



ZW 



The j n is the n-th spherical Bessel function and coefficient can be read from 

n=0 

One easily finds that 
\<Pk(u)\ 2 



sinh(2;'u;) 



(15) 



1 

u 



1 + (mt)' 



' sin 2 (2y) 



2y 2 



0((mty 



(16) 



On the basis of the expansion we see that < <p 2 > is the ultraviolet-finite quantity. It 
is also infrared finite, since the modes are perfectly regular functions for k — > (i.e. 
y->0). 

In the region mt < 1 , with help of the expansion (|T^) , we get the following formula 
for the leading behaviour: 

,2 



< df > 



m 



57T- 



(-ln(mt) + (3/2 - 7 - ln2) + 0((mtf 



(17) 



where 7 is the Euler constant. One should note that the fluctuations vanish as m 
approaches zero and grow with m if we keep mt constant. This agrees with the earlier 
result for an exactly massless field reported in ref. [6]. The interesting feature of 
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the formula flTTD is its non-analyticity in mt and the appearance of the logarithmic 
singularity at t — which is related to the singularity of the RD Universe at t — 0. 

In the regime mt > 1 the numerical study we have performed shows that the integral 
in (0), which is solely a function of mt, has an oscillating behaviour with the period 
close to ir/m and amplitude rising with increasing mt. One can easily check that the 
integral (|H|) is converging rather quickly, hence in a given interval of the variable mt 
one can cut-off the integration from above at a numerically determined value A. Given 
this observation, it is straightforward to find an algebraic approximation for the integral 
as a function of mt. Using the expansion of the Kummer hypergeometric function in 
terms of the modified Bessel functions of the second kind (cf. [7]) 

1, ,1 



1 F 1 (a,b,x) = e2 x T(b - a - -)(-xr " ? 



a-6+i 

2 /v 4 J ' ' 

00 {2b - 2a - !)„(& - 2a) n {b -a-\ + n) 



X £ n\{b) 

X (-l) n Wi+n(-^ 



2 

and expanding the definite integrals in powers of A 2 /mt one obtains in the region 
mt > A 2 the expression of the form 



< (f) 2 > = —y/mt(a sin(2mt) + bcos(2mt) + c + o(A 2 /mt)) (19) 
t z 



where the coefficients a, b, c are 

1 r A * 



a = J^fL dxe ~ 2 ™ MfiWfKx)) (20) 
b = F dxe- 2 ^ (|A(x)| 2 - |/ 2 (a;)| 2 ) (21) 



2(a/2tt) 3 Jo 
1 fA : 



, dxe~ 2 ™ (IfAxM 2 + \f 2 (x)\ 2 ) ^A 3 (22) 

with functions fx, f 2 defined as follows 

r/N tv 1 . ,^(-l/2-2ix) n (l/A-ix + n)(-2ix) n /0 . 1/A . , 
h{x) = T(---ix)J2- T7T, J -cosn/2(n + l/A-tx) 

4 n=0 n \ 2 )n 

_1 ~ (1/2 -2ix)Jl/± - ix + n)(-2ix) n 
+ A^T(--ix)Y J — — -jtt\ — cos7r/2(n + 3/4-ia;) 

4 n=0 n -\2> n 

(23) 
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f2(x) 



i 



oo 



(-1/2 - 2ix) n (l/4 -ix + n)(-2ix), 



n 



sin7r/2(n + 1/4 — ix) 



+ 4 v ^r( i 



4 



-ix) J! 



n=0 



oo 



n!(|)„ 

(1/2 - 2ix) n (l/4 - is + n)(-2ix) 



sin 7r/2(n + 3/4 — is) 



(24) 



where a n = a ra _i(a+n— 1), ao = 1- For the reasonable value of A = 3.5 the least squares 
fit to numerical data in the interval (4.0, 17.0) gives a = —5.9 10~ 4 , b = 4.2 10~ 3 , c = 
1.5 10~ 3 . We note the leading dependence of the < (f 2 > R on m in this range of mt: 
the fluctuations grow proportionally to the square root of the mass. 

One should also note that the fluctuations < <p 2 > R as defined here, see (11), (12), 
are not positive definite. Remarkably, it can easily be shown, that if one perturbs 
the definitions of the functions Ai^{k,m) allowing for the appropriate dependence 
on the ratio which amounts to the modification of our Hilbert space, then the 
regions of an oscillating behaviour with negative values of < cf 2 > R are pushed towards 
increasing values of the argument mt. Since in this letter we do not discuss any specific 
modification of the Hilbert space beyond the simple and most natural definitions (10) 
and (11), the non-positiveness of < <ft 2 > R demands special care when one considers 
physical applications of the present result, as we do in the next Section. The point is 
that a physically meaningful quantity which has the interpretation of the dispersion 
squared should be strictly positive definite. However, we do not rely on the oscillating 
behaviour of < 2 > R in the discussion of applications of the present result being 
only interested in the overall time dependence y/m/t 3 ^ 2 . This leading time and mass 
dependence in the region mt > 1 we believe to be universal, hence we just neglect the 
scheme-dependent oscillating contributions in what follows. 

Up to now we have been calculating curved space vacuum expectation value of 
<p 2 . However, if we were to take into account that the Universe is "hot", i.e. it is 
in fact in a mixed state to which many-particle states may contribute significantly, 
we should better calculate a thermal average of 2 , with finite temperature effects 
included. Assuming thermal equilibrium of the content of the Universe we have 



where <pk are modes given by (|), (|I0"D, and is the occupation number for the particles 




(25) 
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with the comoving momentum k. As we take 

n k = / = = / 1 = (26) 

exp(±Jk 2 /u + m 2 ) - 1 exp(Jk 2 /T 2 + m 2 /T 2 ) - 1 



which is correct for sufficiently large k in view of the choice @. The (|25|) is again 
divergent. However, as usual in finite-temperature calculations, it may be divided 
into T = part and the temperature correction, among which only the former is UV 
divergent. Hence, we can use the renormalization procedure (|TTD to get meaningful 
results even at T > 

< 2 > \renorrnalized,T>0 ~~ > < ^ >_R+ < ^ >JJ (27) 

where 

< 2 >l= 2 / ^!^|0 fc | 2 1 (28) 

J (27r) 3 2A; exp(^k 2 /T 2 +m 2 /T 2 )-l 

This expression may be approximated analytically in two limiting cases: a) ^ >> 1, 
and b) ^ << 1. In the case b) one easily gets 

< 2 >l= ^ (29) 

exactly as in the flat Minkowski case. In the case a) one can see that in the region 
which dominates the integral, k < ^^ a , the second term in ([|) is unimportant. Hence 
we obtain 

< 4> 2 >l= r(5/4)^(^) 3 / 4 -L(T fi /T) 4 ^T 3 e-/ r cos 2 (mt - ?*) (30) 
which is exponentially suppressed. 

2. Implications for particle physics in the expanding Universe 



Let us consider a global U(l) symmetry realized in a single complex scalar field 
model. If Q is the charge of the field Xi the Noether current associated with that 
symmetry is 

f = tQ{xd^x ~ Xd»x} (31) 
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(we put Q=l in what follows) and the conservation law for in the expanding Universe 
reads 

c)V f)V 

d,{a\t)f) = - ia \t){x^ - X^} (32) 



One can see that a symmetry is broken once the rhs of fl32|) is nonvanishing. One can 

see also that when a symmetry is broken explicitly, the net cosmological charge density 

gets generated according to the formula 

, d ,. n ,,,, ., dV dV ' .„„. 
° ^V(t))«-.{X^-X^} (33) 

Let us assume that the term violating the symmetry is 

5V = ^r^(p 2n (34) 

where <f) = Re(x) (in this section we assume the absence of derivative couplings, they 
will be discussed later). Suppose that the initial conditions and the shape of the 
potential are such that the Im(x) and its fluctuations are negligible when compared 
with Re(x) at any time t (this situation may be easily realized in the Affleck-Dine 
model, cf. [8]). Hence 

a- 3 |(A 3 (t))-^0 2 " (35) 
We can see that the magnitude of the symmetry violation is proportional to a coupling 
A, inverse powers of some scale A if n > 2, and to some power of the scalar field 0. 
In this sense one can say that, A and A being fixed in a given theory, it is the <fi what 
determines the amount of symmetry breaking. Here the quantum fluctuations of the 
field <fi come into play. In the quasiclassical picture one can describe the evolution of the 
quantum field, lets call it $, writing it down as the superposition of the quasiclassical 
field cf) which obeys essentially classical (perhaps perturbatively corrected) equation of 
motion and quantum fluctuations 5(f), the dispersion squared of which we identify as 
< 2 > R . If the potential drives the quasiclassical field to zero, then it may happen 
that the magnitude of the symmetry breaking term is determined by the dispersion of 
5<p. Using < > 2n = (< 2 > R ) n one gets an estimate 

«- 3 |(A 3 (t))=^(<0 2 >^r (36) 

Of course, whether this term is really a dominant one or not, it depends on the relative 
magnitude and time-dependence of < (p 2 > R and the classical part of the field. We 
shall investigate this issue later in this work. 
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At this point let us discuss the explicit form of < <fi 2 > R as a function of the tem- 
perature. We know that the time dependence for this quantity is given by ( |TTD if t < — 
and by ( |TT?D if t > — . Let us say assume that the RD epoch starts at t — to. Then 
mt w a 2 (ti) i wnere H(t ) is of the order of the Hubble parameter during inflation, 
H(t ) ~ Hj ~ 10 14 GeV. Assuming then an adiabatic expansion in the RD epoch we 
get 

, 9 m 2 , ,HjT 2 , 

<<P ^-tM-W) ( 37 ) 



R 



as long as T » T R ^- where Tr is the reheating temperature after inflation 

45 4 i i 

T R = (^g-) mffl^Mp)^, (r 7 M P )^] (38) 

(here is the number of relativistic degrees of freedom at of reheating and Tj is the 
total decay width of the inflaton field - a field which drives the transition from the de 
Sitter to RD epoch). Let us take for simplicity the case of a "good" reheating which 
corresponds to Tr « Hj. This gives the condition 

T >% = {mH I y (39) 

If we take m < 10 2 GeV, Hi w 10 14 GeV then we get % < 10 8 GeV. We note that 
for a really soft potential with m ~ 10 -21 eV, one has T* m 10 eV which means that 
in that case the regime where < 2 > R changes only logarithmically extends over the 
whole RD epoch. For T < T* we have instead of 



< 2 >* = ^T 4 (40) 

(g* is the number of relativistic degrees of freedom at temperature T). 

Now, let us investigate the evolution of the classical field (p. For simplicity we 
assume that this evolution is dominated by the mass term in the potential, which is 
usually a good assumption at least in the perturbative regime. In this case the general 
solution to the equation of motion is 

1 



' mt j 4 



T [CJ 1/A (mt) + DJ_ 1/4 (mt)} (41) 



If we set initial conditions at t such that z — mt « 1, then C = r(5/4)0 o + 29t ^° z ° 

m2 3 / 2 



and D = _ 4r ( 3 / 4 )^o z o 1 Q = 0( m t o ) ; Q t( j) = d t (j)(mto)- This gives the <fi at 
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late times, mt > 1, in the form 

« ^(^)~ 3/4 [r(5/4)0 o + 2 -^-\ cos(mt - 3vr/8) (42) 

From previous analysis we have learnt that for T < T* the < <p 2 > R falls off as 

hence we conclude that <( ^,^ R \t<t» ~ m2 2zn&ri , n - which does not depend on time 

4> 2 i 1 ^ 1 * [r(5/4)<^o+ °^ tTO ] 2 

(and temperature). This implies that if the initial conditions happen to make this ratio 
large, then the fluctuations dominate over the classical part of the field. Moreover, in 
general the field <fi has some additional couplings to light particles, which facilitate decay 
of the field with the decay width T^. This changes the behaviour of the classical field 
(J), namely <fi 2 — > exp(— Ft) cf 2 . Actually, as pointed out by several authors in the context 
of the Affleck-Dine mechanism (which corresponds to our toy model when n=2) the F^ 
should be large in order to avoid an unobservable excess of the net charge produced 
during symmetry violation [8]. We want to stress that in such a case, the < 4> 2 > R , 
decaying accordingly to the power law, dominates the divergence of the Nother current 
and the net cosmological charge density even at the late times. 

Let us check whether at T > T* the < <fi 2 > R may be significantly large. From (57) 
and ( [12]) we obtain, averaging over oscillations, 

< d> 2 > R , 2m 2 



It>t *~ (r(5/4)0 o + ^Mo; 2 



(43) 



H < m (44) 



That means that fluctuations are important as long as the condition 

z d t (j)o 
m 

is fulfilled. One can see that even if m « Hi and 0o ~ Hi, which happens to be 
the case if the initial conditions at t — to are produced by large quantum fluctuations 
in the preceding de Sitter epoch [8], the condition ( fP|) may be fulfilled provided that 
d t (po is sufficiently large and negative. 

Finally, let us consider models where a massive scalar <ft is derivatively coupled to 
other particle species. This situation corresponds for instance to models possessing 
pseudogoldstone bosons with nonvanishing masses. The relevant scenario is similar 
to that of the "spontaneous baryogenesis" described in ref. [5]. If a Lagrangian has a 
coupling of the form = —j4>d^j^ (f being some, presumably large, mass scale) where 
d^j^ is a divergence of a current corresponding to some explicitly broken symmetry, the 
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baryon number symmetry for instance. Then, as we have shown, there are fluctuations 



in the field with dispersion y < 2 > R . We may represent them as the effective term 
in the Lagrangian 

1 



f v < 2 > R d lx f (45) 



Up to the total divergence fl4"5| ) is equivalent to 

1 



d^<^ >B j G (46) 



This produces an effective chemical potential \x = — j<9oy< 4> 2 > R for the charge den- 
sity j° which means a nonzero cosmological charge density generated in thermal equi- 
librium. Explicitly, cf. [9], 



1 



3° « -jdty/< 2 > R T 2 (47) 



or charge to entropy ratio 



f/s K -—dty/< 2 > fi (46 



where is the number of relativistic degrees of freedom at temperature T. The above 
estimate gives in the case of our toy model 

o m 1 

3 /S * 49*fT ty /]n(l/mt) 1 ' 

for T > and 

for T < T*. One can see that both expressions fall off as time elapses, and that the 
decrease at T < T* is faster than at T > T*, essentially ~ T 5//2 below and ~ T 
above. If there is no phase transition in the model before the end of RD epoch, then 
the final charge to entropy ratio produced will be equal to (|49|) or (]50|) taken at the 
"decoupling" temperature Tjj. This is the temperature at which symmetry violating 
interactions fall off from equilibrium or the one which corresponds to the end of RD 
stage, when the shape of the fluctuations changes qualitatively i.e. at To rs Tf close 
to 10 eV. As previously, the numerical values predicted depend on various details of 
a model under investigation. For example, let us take Tp = 10 eV and = 100. 
Then if we require the charge-to-entropy ratio to be equal to 10 -10 , as it should be 
for the baryonic charge, then we get the condition m = f A x 10~ 77 GeV, which gives 
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m = 10" 17 GeV for / = 10 15 GeV and m = lGeV for / = 10 19 GeV. 



3. Conclusions 

In this letter we have found explicit expressions for a massive scalar field fluctua- 
tions in the flat radiation dominated universe. It turns out that in the region of small 
mt, i.e. shortly after the beginning of the RD epoch or for very light fields, the fluc- 
tuations decrease with time only logarithmically and are proportional to the square 
of the mass of the field in question. For large mt, i.e. very late or for a heavy field, 
the time dependence is stronger, l/t 3//2 , but the mass dependence becomes weaker - 
proportional to the square root of m. As far as finite temperatures are concerned, 
we have concluded that the "radiation-dominated" background modifies Minkowski 
space results rather weakly. At low temperatures, i.e. at large ratios m/T, the ther- 
mal contribution is exponentially suppressed the suppression becoming stronger as the 
temperature decreases with time. At high temperatures the result coincides essentially 
with that of Minkowski space. In general, fluctuations vanish when one takes the limit 
m — > 0. 

Given all that we argue that the fluctuations may still play a significant role in 
particle physics models, which has been illustrated in the second part of the work. 
Within the family of models we discuss, the case when our parameter n equals 2 
corresponds precisely to the Affleck-Dine model, and the higher n terms are often 
encountered in the important class of string inspired models. Hence we conclude that 
fluctuations we have described constitute the phenomenon which is relevant in a very 
general situation when some cosmological charge density, first of all the baryonic charge 
density, is supposed to be generated during the radiation dominated epoch. We also 
note that although the inflationary scenario is widely accepted, our results do not rely 
on the existence of the de Sitter epoch preceding the RD stage in the early Universe. 

In conclusion, we have demonstrated that quantum fluctuations due to a nontriv- 
ial gravitational background during radiation dominated epoch in the evolution of the 
Universe do in fact exist and may have observable consequences for cosmology of real- 
istic particle models. 
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